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Pumping velocity in homogeneous helical turbulence with shear 
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Using different analytical methods (the quasi-linear approach, the path-integral technique and tau-relaxation 
approximation) we develop a comprehensive mean-field theory for a pumping effect of the mean magnetic field 
in homogeneous non-rotating helical turbulence with imposed large-scale shear. The effective pumping velocity 
is proportional to the product of a effect and large-scale vorticity associated with the shear, and causes a separa- 
tion of the toroidal and poloidal components of the mean magnetic field along the direction of the mean vorticity. 
We also perform direct numerical simulations of sheared turbulence in different ranges of hydrodynamic and 
magnetic Reynolds numbers and use a kinematic test-field method to determine the effective pumping velocity. 
The results of the numerical simulations are in agreement with the theoretical predictions. 



PACS numbers: 47.65.Md 



I. INTRODUCTION 



The origin of cosmic magnetic fields is one of the funda- 
mental problems in theoretical physics and astrophysics. It is 
generally believed that solar and galactic magnetic fields are 
caused by the combined action of helical turbulent motions 
of fluid and differential rotation |[ll-|7|- In most of these stud- 
ies, differential rotation plays merely the role of enhancing the 
magnetic field in the toroidal direction. However, in recent 
years there has been increased interest in mean-field effects 
caused specifically by turbulent shear flows. This interest is 
caused by discoveries of the shear dynamo [8,9] and vortic- 
ity dynamo lliol flld in non-helical homogeneous turbulence 
with a lar ge-s cale shear. In particular, recent numerical ex- 
periments lfl2l - [l7ll have clearly demonstrated the existence of 
a shear dynamo of a large-scale magnetic field in non-helical 
turbulence or turbulent convection with superimposed large- 
scale shear However, the origin of the shear dynamo is still 
subject of active discussions ifsllgl [Tsi [T8I - I23I1 . 

There are three additional phenomena that are also related 
to the presence of shear One is the vorticity dynamo, which is 
the self-excitation of large-scale vorticity in a turbulence with 
large-scale shear It has been predicted theoretic ally lUOl [Till 
and detected in recent numerical experiments 1IT3I [T4l 12411 . 
The vorticity dynamo can also affect the dynamo process of 
the mean magnetic field. Another phenomenon is a non-zero 
a effect in non-helical turbulence with shear when the system 
is inhomogeneous or density stratified. In that case there is 
an a effect ^ [Tsll that can lead to an alpha-shear dynamo. 
Finally, when homogeneous turbulence with shear is heli- 
cal, there is an effective pumping velocity 7 oc aW of the 
large-scale magnetic field, where W is the large-scale vortic- 
ity caused by shear This effect has so far only been found 
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in direct numerical simulations (DNS) |[25ll . but there has so 
far been no theory for this new effect, nor has there been a 
systematic survey of DNS for determining the dependence of 
pumping on magnetic Reynolds and Prandtl numbers as well 
as the turbulent Mach number 

The goal of the present study is to develop a comprehensive 
theory of mean-field pumping in homogeneous helical turbu- 
lence with shear and to perform systematic numerical simu- 
lations designed for detailed comparison with the theoretical 
predictions. It is important to emphasize that the pumping of 
the large-scale magnetic field discussed usually in the litera- 
ture has always been connected with inhomogeneous turbu- 
lence iHH^Il^j, but here we study the pumping for homoge- 
neous, albeit helical turbulence. 

II. GOVERNING EQUATIONS 

We consider homogeneous helical turbulence with a linear 
shear velocity U = (0, Sx, 0). Averaging the induction equa- 
tion over an ensemble of turbulent velocity field yields the 
mean-field equation: 



dB 



-— = Vx (UxB + uxb-riVxB) , 

at ' 



(1) 



where Ei = {u x b)i = aijBj + bijkVkBj is the mean 



electromotive force, u and b are the fluctuations of velocity 
and magnetic field, overbars denote averaging over an ensem- 
ble of turbulent velocity fields, B is the mean magnetic field, 
U is the mean velocity that includes only the imposed large- 
scale shear, and 77 is the magnetic diffusion due to electrical 
conductivity of the fluid. Note that the part aijBj in the ex- 
pression for the mean electromotive force determines the ef- 
fective pumping velocity, 7^ = —-^eijkCiij, and the a tensor, 

aij = ^{aij + aji), i.e., s'f''' = OLijBj + (7 x B)i, while 
the turbulent magnetic diffusion and the shear-current dynamo 
effect are associated with the bijk term. 

To determine the turbulent transport coefficients in homo- 
geneous helical turbulence with mean velocity shear we use 



2 



the following equations for fluctuations of velocity and mag- 
netic field: 



O 11 

-{U-V)u-{u-V)U--Vp+—\ib-V)B 

p Airp 



dt 



db 



(B-V)u - {u-V)B + (6- V)C7 - (U-V)b 



-r]Ab + b 



N 



(2) 



(3) 



where ly is the kinematic viscosity, p is the mean density of the 
incompressible fluid flow, p is the fluctuation of total (hydro- 
dynamic and magnetic) pressure, the magnetic permeability 
of the fluid is included in the definition of the magnetic field, 
and b^ are the nonlinear terms, and p/^"^ is the stirring 
force for the background velocity fluctuations. 

We begin by deriving expressions for the pumping effect 
that are valid in different regimes, where fluid and magnetic 
Reynolds numbers are both small, both are large, or only the 
fluid Reynolds number is large, but the magnetic Reynolds 
number is small. These results will then be compared with 
those of DNS in the corresponding regimes. 



A. Small magnetic and hydrodynamic Reynolds numbers 

We use the quasi-linear or second order correlation ap- 
proximation (SOCA) applied to shear flow turbulence (see 
III8II20I1 ). This approach is valid for small magnetic and hy- 
drodynamic Reynolds numbers. To exclude the pressure term 
from the equation of motion (|2]i we calculate Vx(Vxm), 
then we rewrite the obtained equation and Eq. (|3) in Fourier 
space, apply the two-scale approach (i.e., we use large- 
scale and small-scale variables), neglect nonlinear terms in 
Eqs. (|2]i-(l3]i, but retain molecular dissipative terms in these 
equations. We seek a solution for fluctuations of velocity and 
magnetic fields as an expansion for weak velocity shear: 



b = 6(0) + 



(4) 
(5) 



where 



dkn 



(6) 



-"(fc, = G^{k,uj) i2kiqSjp + Sij kq-gj^ 

x{VpUq)uf(k,L0), 



(7) 



uf\k,Lu) 



6-k — 



(8) 



Here Gi/(fc,a;) = (vk"^ — iuj)~^ , G,,(fc,aj) = (rjk'^ — iuj)~^ , 
and Sij is the Kronecker tensor The statistical properties 
of the background velocity fluctuations with a zero large- 
scale shear, u^^\ are assumed to be given. For derivation of 
Eqs. ©-(O we use the identity 



Uq{Q) 6„(fc -Q)dQ^ i{VpUq) 



dbn 
dkn 



which is valid in the framework of the mean-field approach, 
i.e., it is assumed that there is scale separation. Equations (|6]l- 
(|8]l coincide with those derived by [18], and they allow us 
to determine the cross-helicity tensor g 



(1) 



,(1) 



bf^). This procedure yields the contributions £m^ 



£mij J gij\k,uj)dkduj to the mean electromotive force 
caused by sheared helical turbulence. We are interested first 
of all in the contributions to the mean electromotive force 
which are proportional to the mean magnetic field, i.e., S^""^ = 
ctijBj + (7 X B)i. For the integration in w-space and in k- 
space we have to specify a model for the background shear- 
free hehcal turbulence (with B = 0), which is determined by 
equation: 



{ui{k,uj)uj{-k,~LL>)) 



,2\(0) 



fc2 



87rfc2 
ki (tt- (V X tt)>(0) 



~fc2" 



(9) 



where E{k) is the energy spectrum (e.g., a power-law spec- 
trum, E{k) cx (fc/fct)~* with the exponent 1 < q < 3 
for the wavenumbers kf < k < fcj, where kf and k^ are 
the forcing and dissipation wavenumbers), and Sijk is the 
fully antisymmetric Levi-Civita tensor We consider the fre- 
quency function ^{oj) in the form of the Lorentz profile: 
^{lj) ~ vk!^ /[n (w^ + v^k'^)]. This model for the frequency 
function corresponds to the correlation function 

(ui(t)uj(i + r)) cx exp(-r t/fc^). (10) 

In that case, and under the assumption of small magnetic and 
hydrodynamic Reynolds numbers, the effective pumping ve- 
locity, 7, and the off-diagonal components of the tensor cky 
are given by 



7 = 



Ci{q) 



Pm 



Oil 



2 V 1 + Pm 
Ci(q) (2Pm + l)Pm 



Re^ Tf a* W, 



(11) 



Cl{q) 



(1 



Pm)2 

-4 



ReVfa, (5C/)„-, (12) 



dk 



bf\k,uj) = G,,(fc, I i{k-B)Sij ~ (s.^j km-^ where a 



q-1 



k 
ki 

1 - {kf/k^Y+^ 



1 - ihlk^Y' 



(13) 



-(Tt/3) (tt • (V X u)Y^\ Pm = v/r] is the 
magnetic Prandtl number. Re — Tf {u'^)^'^^ /v is the hydro- 
dynamic Reynolds number, Rm = RePm is the magnetic 
Reynolds number, Tf = If /wrms is the turnover time, where 
If = 1/fcf is the energy-containing (forcing) scale of a ran- 
dom velocity field, and Wj-ms = \/(u^^- For the integration 
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in w-space we use the integrals given in Appendix lAl 

For linear sheai" velocity, U — (0, Sx, 0), the mean vorticity 
is = V X t/ = (0, 0, S), and the mean symmetric tensor 
{dU)ij = iU j + V jUi)/2 has only two nonzero compo- 
nents: {dU)i2 = {dU)2i — S/2. Therefore, has two 
non-zero off-diagonal components caused by both, shear and 
helical turbulence ai2 = a2i, while the effective pumping ve- 
locity, 7, has only one component directed along the vertical 
axis, 7 ~ (0, 0, 7): 



Ci(g) 



Pm 
1 + Pm 
Cijq) (2Pm 
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Re^ a* 



Sh, 



(1 + Pm) 



l)Pni 2 

He a. 



(14) 



Sh, (15) 



where Sh — rt S is the shear parameter As follows from 
Eqs. ([Til and 7 cx Pm and cx Pm for Pm <^ 1, 
while for Pm ^ 1 the effective pumping velocity 7 and 
are independent of Pm. For all values of the magnetic Prandtl 
numbers, 7 and ai2 are positive. This asymptotic behavior 
which is valid for Re ^ 1, is in agreement with Figs.[T]and|2] 
(see Sect. III). Note that the diagonal components of the ten- 



sor aij in this case are 



C2(g) 



C2iq) 



Rm 
1 + Pm 



Tf (it • (V X u)) 



(0) 



dk 



1 - (fcf/fcd 



1 - (ki/kd) 



q-l 



(16) 



(17) 



B. Large magnetic and hydrodynamic Reynolds numbers 

To determine the the effective pumping velocity and the 
tensor in homogeneous helical turbulence with mean ve- 
locity shear for large magnetic and hydrodynamic Reynolds 
numbers we use the procedure which is similar to that applied 
in in earlier investigations of shear flow turbulence. Let 
us derive equations for the second moments. We apply the 
two-scale approach, e.g., we use large scale R = {x + y)/2, 
K = ki + k2 and small scale r = x — y, k = (fci — fc2)/2 
variables (see, e.g., |28]). We derive equations for the follow- 
ing correlation functions: 



L{ui;uj), hij{k) = L{bf,bj), 



where 



L(a; c) 



[a{k + K/2)c{-k + K /2)) exp {iK-R) dK, 



and (...) denotes averaging over ensemble of turbulent veloc- 
ity field. The equations for these correlation functions are 
given by (see i^) 

dfijik) 



dt 



dhij{k) 



|P = -z(fc.B)<l>„+4 



E. 



S (Tt- 



{U)K 



dt 



iik.B)[f,,ik) - h,,{k) - h\P] 

+ JijmniU)9mn+Sfgij, 



(18) 



where hereafter we omit the arguments t and R in the correla- 
tion functions and neglect small terms ^ O(V^). Here Eij is 
related to the forcing term and V = d/dR. In Eqs. ( fTSI ). 
<^^j{k) = {Ai:-p)~^[g,j{k) - gj,{-k)l and TVA/ij, Nhij, 
Mgij , are the third-order moments appearing due to the non- 
linear terms which include also molecular dissipation terms. 
The tensors lfj„,„iU), (C7) and Jfj,„,,(U) are given 

by 



-^ijmni^) ~ y^irn^jq^pn Sjm^iqSpn 

d \ - 



'^ijmni^) — ( ^kjq^im^pn ^im^pn^ jq ~^ ^jji^pm^ii 



d 



where kij — kikj/k^. The source terms I^^, and If^ which 
contain the large-scale spatial derivatives of the magnetic field 
B, are given in @]. Next, in Eqs. ( fTsT i we split the tensor for 
magnetic fluctuations into nonhelical, hij , and helical, h\^^ , 
parts. The helical part of the tensor of magnetic fluctuations 
h\^^ depends on the magnetic helicity and it follows from 
magnetic helicity conservation arguments (see, e.g., ll29l - l32ll 
and [7] for a review). 

The second-moment equations include the first-order spa- 
tial differential operators J\f applied to the third-order mo- 
ments A/ (I"). A problem arises how to close the system, 
i.e., how to express the set of the third-order terms MM'^™-^ 
through the lower moments M^^^^ (see, e.g., fss'-^s']). We use 
the spectral r-closure-approximation which postulates that the 
deviations of the third-moment terms, MM^^^^^ (fc), from the 
contributions to these terms due to the background turbulence, 
_/\/'7\/(iii.o) ^^^^ expressed through the similar deviations of 
the second moments, Af(")(fc) - M("^°)(fc): 

UM^™\k)-MM^™'°\k) 

= [M(")(fc)-M("-")(fc) 

Tr(k) L 

(see, e.g., ll33l l36i IjtII). where Tr{k) is the scale-dependent 
relaxation time, which can be identified with the correlation 
time, Tf, of the turbulent velocity field for large hydrodynamic 
and magnetic Reynolds numbers. The quantities with the su- 
perscript (0) correspond to the background shear-free turbu- 
lence with a zero mean magnetic field. We apply the spectral 



(19) 
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T approximation only for the nonhelical part hij of the tensor 
of magnetic fluctuations. Note that a justification of the r ap- 
proximation for different situations has been performed in a 
number of numerical simulations and analytical studies (see, 
e.g., Iii8,-45J). 

We take into account that the characteristic time of variation 
of the magnetic field B is substantially longer than the corre- 
lation time Tf . This allows us to obtain a stationary solution 
for Eqs. (fTsT l for the second-order moments, M(")(fc), which 
are the sums of contributions caused by shear-free and sheared 
turbulence. The contributions to the mean electromotive force 
caused by a shear-free turbulence and sheared non-helical tur- 
bulence are given in |9]. In particular, the contributions to 
the electromotive force caused by the sheared turbulence read: 

— £mji J glp (fc) dk, where the corresponding contri- 

( s) 

butions to the cross-helicity tensor g^j in the kinematic ap- 
proximation, are given by 



J:. 



ijmn 



(fc) {k.B) 



-Tr{k) {k-B)I, 



s 

ijmn 



f(0) 



(20) 



and we use the following model for the background shear-free 
helical turbulence (with B — 0): 



/(°) = K(fc)t.,(-fc,)>(") 



ki (tf(VXM))(0) 



m 

87rfc2' 



(21) 



where the energy spectrum is E{k) — {q — I) (k/kf)^'^, 
k[ = l/£f and the length £{ is the maximum scale of tur- 
bulent motions. The turbulent correlation time is Tr{k) = 



2Tf (fc/fcf 



Therefore, for large magnetic and hydrody- 



namic Reynolds number the effective pumping velocity, 7, 
and the off-diagonal components of the tensor aij caused by 
sheared helical turbulence are given by 



2 



^-(5~2g)Tfa, (dU), 



(22) 
(23) 



Since the mean symmetric tensor {dU)ij has only two 
nonzero components: {dU)i2 — {dU)2i — S/2, the tensor 
aij has only two non-zero off-diagonal components, 0:12 = 
a2i- In particular. 



7 = - a* Sh, 



(24) 



2 2 
ai2 = "21 = - - (5 - 2q) a, Sh -- a* Sh, (25) 
5 

where we have used the Kolmogorov kinetic energy spectrum 
exponent q ~ 5/3 in Eq. ( |25] |. The diagonal components of 
the tensor in this case are a — a^, (see, e.g., 1 1, 3]). These 
results for large magnetic and hydrodynamic Reynolds num- 
ber are in qualitative agreement with DNS performed in [|25l . 



C. Large magnetic Reynolds numbers and small 
hydrodynamic Reynolds numbers 

To develop a mean-field theory for large magnetic Reynolds 
numbers and small hydrodynamic Reynolds numbers we use 
stochastic calculus for a random velocity field. To derive an 
equation for the mean magnetic field we use an exact solu- 
tion of the induction equation for the total field B (which is 
the sum of the mean B and fluctuating b parts) with an initial 
condition B{t — to,x) = B{to,x) in the form of a func- 
tional integral: 

B,{t,x) = {G,j{t,to,0 exp{i-V)Bj{to,x))^, (26) 



(see, e.g., 1461 147|]), where the operator exp(4 . V) is deter- 
mined by 



00 

exp(f V) = 5: - (I . 



k=0 



(27) 



^ = ^ — x (see AppendixlBt. The Wiener trajectory ^{t, to, x) 
is determined by 

ta,x) = x- f v{U, I) da + {2rif/^w{t - to), 
Jo 

(28) 

where t^r = t — a, and the velocity field v is the sum of 
the mean shear velocity U and fluctuating u parts. We con- 
sider large magnetic Reynolds number, but take into account 
small yet finite magnetic diffusion 77. The magnetic diffu- 
sion can be described by a random Wiener process w{t) that 
is defined by the following properties: {wi{t))w — and 
{wi{t + T)wj {t))w ~ ^Sij, where {■)^ denotes the averaging 
over the statistics of the Wiener random process. The function 
Gij [t, s, ^) is determined by equation: 

dG,,it,s,^) 



ds 



N,kGkj{t,s,^), 



(29) 



Sij and Nij — VjVi. 



with the initial condition Gij (t = i 

The form of the exact solution (|26] l allows us to separate the 
averaging over random Brownian motion of particles (i.e., the 
averaging over a random Wiener process w{t)) and a random 
velocity u. 

We consider a random flow with a small yet finite Strouhal 
number (that is the ratio the correlation time of a random fluid 
flow to the turnover time If/umis)- A random velocity field 
with a small Strouhal number can be modelled by a random 
velocity field with a constant renewal time t. Assume that in 
the intervals . . . (— r, 0); (0, r); (r, 2t); . . . the velocity fields 
are statistically independent and have the same statistics. This 
implies that the velocity field looses memory at the prescribed 
instants t = mr, where to = 0, ±1, ±2, . . .. This velocity 
field cannot be considered as a stationary velocity field for 
small times ^ r, however, it behaves like a stationary field for 
t ^ T. Averaging Eq. (|26l l over the random velocity field we 
arrive at the equation for the mean magnetic field, B: 

OB ■ 

—l=[yx (UxB)]^ + A,,^VmBj 



5 



(see AppendixlBl). where 



T 



D. 



IT 



(31) 
(32) 



the angular brackets (•) denote an ensemble average over the 
random velocity field. Therefore, the mean magnetic field is 
determined by double averaging over two independent ran- 
dom processes, i.e., by the ensemble average over the random 
velocity field and by the average over Wiener random process 
w{t). 

We are interested in the lowest-order contributions to the 
mean electromotive force which are proportional to the mean 
magnetic field, ff"-* = a^.J Bj, where a^j = (l/2)£i„„ Anjm 
and the tensor Aij„i reads: 



4 ■ ■ — — 



i dt £ dt' ( i)] ^ [v,v,{t', i)]^), 



(33) 



where x ^ y and [w,„(t,^)]^ denotes the Eulerian velocity 
determined at the Wiener trajectory ^ that passes through the 
point X at instant t. Hereafter the angular brackets denote 
double averaging over a random velocity field and over the 
statistics of the Wiener process. 

For small hydrodynamic Reynolds numbers we seek the so- 
lutions of the linearized Navier-Stokes equation Q for incom- 
pressible velocity field u as superpositions of the Orr-Kelvin 
random shearing waves u{t,r) — J u{t,ko) exp[ik{t) ■ 
r]dko, where feo = (k j,p,k,,,K ), k{t) = {kj^o - 
Skyt,ky,k^) (see, e.g., BMjSCT). Therefore, the effec- 
tive pumping velocity, 7, and the off-diagonal components of 
the tensor ay are given by 



1 _ 1 



Z 

4^ JO 

X kmit') (u™(t,fcoX(i',fco)), 
I , .1 



dt 



dt' 



(34) 



-jnm ^nim ) 



dt (Sinm kj {t ) -'r 



x{v„i{t,ko)vl{t' ,ko)). 



jnm ki(t ) ) 

(35) 



Using these equations and Eqs. (|C13t - (IC10l l in Appendix ICl 
we obtain the effective pumping velocity, 7 = (0, 0, 7), and 
the off-diagonal components ai2 — a2i of the tensor aij for 
large magnetic Reynolds numbers and small hydrodynamic 
Reynolds numbers: 



■ a, Sn Re , 



ai2 



3 

"21 



(36) 



C2(g)Re- 



Tf 



a* Sh, 



(37) 



where Re ^ t/t{ < 1. The diagonal components of the 
tensor in this case obtained using path-integral approach 



are a = -(1/3) {tu ■ (V X (see, e.g., iBilll). In the 

next section we discuss comparison with new systematic DNS 
designed for comparison with our theoretical predictions. 



III. COMPARISON WITH DNS 

A. Numerical model 

Our DNS model is identical to that used in ||25|] . We be- 
gin by testing the analytical results numerically using three- 
dimensional simulations of isotropically forced turbulence 
in a fully periodic cube of size {^tt)^. The uniform shear 
U = (0, 5*2;, 0) is imposed using the shearing box method and 
the gas obeys an isothermal equation of state characterized by 
the constant speed of sound Cg. We solve the continuity and 
Navier-Stokes equations in the form 



Vlnp 
Vt 



= -U-V\np-V-U, 



(38) 



VU 

'vt 



= -U-VU -SU^y-clV\np + f + F^,,„ (39) 



where the imposed shear is subsumed in the advective deriva- 
tive 



V _ d d 
Vt dt dy 



(40) 



Here p is the density, U is the velocity, / describes the forc- 
ing, and i^visc = /o~^ V • {2pvS) is the viscous force, where v 
is the kinematic viscosity, and 



(41) 



is the traceless rate of strain tensor The forcing function / is 
given in [52]: 

f{x, t) = Re{Nfkit) exp[ifc(0 • x + i0(i)]}, (42) 

where x is the position vector The wavevector k{t) and the 
random phase — tt < (f){t) < tt change at every time step, so 
f{x,t) is (5-correlated in time. The normalization factor N is 
chosen on dimensional grounds tobe = fQCs{\k\cs/Sty/'^, 
where /□ is a nondimensional forcing amplitude. At each 
timestep we select randomly one of many possible wavevec- 
tors in a certain range around a given forcing wavenumber 
The average wavenumber is referred to as kf. In the present 
study we always use kf/ki = 5. We force the system with 
transverse helical waves 05311. 



/fc = R-/i"°'^'^ with R^^_^!^_J^!M|^, (43) 

where a = 1 for the fully helical case with positive helicity of 
the forcing function. 



/(nohcl) 
k 



(fc X e) /v/fc2 - (fc . e)2 



(44) 



6 



is a non-helical forcing function, and e is an arbitrary unit 
vector not aligned with k; note that — 1. We use fully 

helical forcing, i.e. cr = 1, in all of our runs. 

The boundary conditions in the y and z directions are peri- 
odic, whereas shearing-periodic conditions are used in the x 
direction. The simulations are governed by the fluid and mag- 
netic Reynolds numbers, the magnetic Prandtl number, and 
the shear and Mach numbers: 



Re 
Sh 



vkf 
S 



Rm = 
Ma 



■qkf 

^rms 



Pm 



(45) 



Here Wj-ms is the root mean square velocity of turbulent mo- 
tions and rj is the magnetic diffusivity. We use the PENCIL 
to perform the simulations. 



a 



10" 



10" 



10" 



10" 



10" 



10" 



0.01 



Re = 0.16, C,(g) = 1.05 ; 



' Re = 0.04, C,(qr) = 1.2 



0.10 



1.00 
Pm 



10.00 



B. Test field method 

We apply the kinematic test-field method (see, e.g., ifTsl 
I54. 55,1) to compute the effective pumping velocity, 7, and all 
components of the tensor aij . The essence of this method is 
that a set of prescribed test fields B'^p^''^ and the flow from the 
DNS are used to evolve separate realizations of small-scale 
fields b^P'i\ Neither the test fields B^^ "?) nor the smafl-scale 
fields act back on the flow. These small-scale fields 

are then used to compute the electromotive force cor- 
responding to the test field The number and form of 
the test fields used depends on the problem at hand. For the 
purposes of the present study we use uniform horizontal test 
fields 5(1) = (So, 0,0) and!B(2) = (0, Bq, 0), in which case 
the series expansion of the electromotive force contains only 
a single term 

Ef'^a^^B,. (46) 

We present the results using the quantities: 

a = 5(011 + 022), (47) 
ai2 = a2i = 5(021 + 012), (48) 
7= 5(021 -012). (49) 

We use ao — as a normalization factor when present- 

ing numerical results. Errors are estimated by dividing the 
time series into three equally long parts and computing time 
averages for each of them. The largest departure from the 
time average computed over the entire time series represents 
the error This definition of the error bar gives an indication 
about the mean value that one would obtain for shorter parts 
of the time series. With this definition, the error bars do nor- 
mally become shorter for longer runs, provided the time seri es 
is stationary. This would not be the case for the rms value of 
the deviations, which might sometimes also be of interest. 
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FIG. 1. Pumping coefficient 7 = ^(ffl2i — Q.12) normalized by qq = 
iurms as a functioH of Pm for two values of Re (Sets Al and A2). 
The shear parameter Sh = -0.20 (-0.13) for Re = 0.04 (0.16). 
Analytical results according to Eq. ( I14t are overplotted with dotted 
lines. The values of Ci (g) are used as fit parameters and indicated in 
the legends. 



C. Results 

We perform several sets of simulations where we vary the 
parameters Pm, Rm, Sh, and Ma individually to study the 
analytical results derived in Sect. II; see Table J] The setup 
used here is prone to exhibit the so-called voiticity dynamo 
ifioi I24I1 . due to which large-scale voiticity is generated, and 
complications can arise in the interpretation of the simulation 
data. Here we restrict the studied parameter range so that the 
values of Re and Sh are subcritical for the vorticity dynamo. 
In our runs where the Reynolds numbers are of the order of 
unity or less, a low grid resolution of 32'^ is often sufficient. 
Indeed, in Table II we show the results obtained for differ- 
ent resolutions ranging from 16"^ to 128"^ for Rm around 1, 
which demonstrates good convergence of the results within 
error bars. 



1. Dependence on Pm 

Figure [U shows our results for 7 as a function of magnetic 
Prandtl number Pm. We find that the numerical results coin- 
cide with the analytical formula, Eq. ( fT4l i. Values of the order 
of Ci (q) « 1 fit the DNS results within the error estimates. 

Figure |2] shows the results for 0112 as a function of Pm for 
two values of Re. The data for Q!i2 shows significantly larger 
fluctuations than the corresponding results for 7. However, the 
DNS results seem to fall in line with the analytical expression, 
Eq. ( fTSl ). although the value of Ci (q) needed to fit the data is 
an order of magnitude larger than in the case of 7. This can 
be explained by comparing Eqs. (|36] | and (IJTI i. which show 
that 7 cx Re^, while q;i2 oc Re(r/rf), where r is the flow 
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TABLE I. Summary of the runs. 



t>et 


Ke 


Pm 




on 




Ma 


grid 


A 1 

Al 


U.U4 


U.U5. . . 


23 


n on 
— U.2U 




n m n 


oo3 c /1 3 




U.lO 


U.Uz. . . 


on 


n 1 Q 




n m 
U.UiD 


61 . . . d4 




n no o 1 






n no K 




n non 


oo3 OKt;3 


B2 


0.08. . . 83 






—0.075 




0.080 


o no o r" /^3 

32 . . . 256 


B3 


0.08. ..3.5 






-0.25 




0.080 


32^ 


B4 


0.08... 0.4 






-2.5 




0.080 


32^ 


CI 


0.04 






-0.020... - 


0.19 


0.010 


32^ 


C2 


0.16 






-0.012... - 


0.12 


0.016 


32^ 


C3 


0.45 






-0.009... - 


0.09 


0.023 


32^ 


C4 


1.3 






-0.006... - 


0.07 


0.032 


32^ 


Dl 


0.08 






-0.010 




0.002 . . .0.41 


32^ 



TABLE n. Convergence study of 7 and Q21 for Rm = 1.3 and 
Sh = —0.06 from simulations with different grid sizes. 

Run 7/ao[10~^] 0:21/0:0 [10~^] grid 



El 1.02 ±0.12 0.94 ±0.25 16'' 

E2 1.05 ±0.07 0.89 ±0.20 32^ 

E3 0.99 ±0.06 0.83 ±0.55 64^ 

E4 0.94 ±0.18 0.88 ±0.23 128^ 



10" 



10" 



10" 



10" 



Re = 0.16, C,{g) = 8 



Re = 0.04, C,(g) = 22 



0.01 



0.10 



1.00 



10.00 



Pm 



FIG. 2. Symmetric part of a^j, 012 = ^(021 ± ai2) normalized by 
Oo = fwrms as a function of Pm for the same runs as in Fig.[T] The 
dotted hnes show the analytical result according to Eq. ( I15t . with the 
values of Ci (g) indicated in the legends. 



renovating time, and Tf = ^f/wrms is the turnover time of 
turbulent eddies. Note that Eqs. (|36] | and (l37l l are obtained for 
large magnetic Reynolds numbers, while Re ^ r/rf < 1. 
This implies that for these conditions 0:12 ^ 7. The latter is 



0.100 



f 0.010 



0.001 



Re = 0.16, Cj^q) = 1 



Re = 0.04, Cj(g) = 0.15 



0.01 



0.10 



1.00 



10.00 



Pm 



FIG. 3. o-effect as a function of Pm normalized by 00 — 3 Unns for 
the same runs as in Fig.[T] Analytical results according to Eq. l il6t 
are overplotted with dotted lines. The values of C2{q) are used as fit 
parameters and indicated in the legends. 



in agreement with DNS results (see Figs.[T]and|2]l. 

In Fig. [3] we show a-effect (the diagonal elements In Fig. [3] 
we show the a-effect (the diagonal elements of the aij tensor) 
as a function of the magnetic Prandtl number Pm. These re- 
sults are in a good agreement with the analytical results (fTSI ). 



2. Dependence on Rm 

Our results for 7 as a function of Rm are shown in Fig. |4] 
We find that for Rm smaller than roughly two, 7 is well de- 
scribed by the analytical result, Eq. (fl4] i obtained for Rm <C 1 
and Re ^ 1. For greater Rm, 7 is consistent with a constant 
value as a function of Rm, and is in accordance with Eq. (|24] | 
derived for Rm ^ 1 and Re ^ 1. Note also that for the largest 
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0.1 



1.0 10.0 
Rm 



100.0 



1.0 



0.1 



Sh = -0.075 
Sh = -0.025 



// 
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/ 



1.1 

I 

^ 0.8 



0.1 1.0 10.0 100.0 
Rm 



0.1 



1.0 10.0 
Rm 



100.0 



FIG. 4. 7 as a function of Rm for Pm = 1 and for four values of 
Sh (-0.025, -0.075, -0.25, and -2.5; see Sets Bl to B4). The 
hnes show the analytical results according to Eqs. ( 114) and l l24t with 
Cx{q) = 1, for Sets Bl (dotted hnes), B2 (dashed), B3 (dot-dashed), 
and B4 (triple-dot dashed), respectively. 



FIG. 6. a-effect as a function of Rm normalized by qq = ^Mrms for 
two values of Re (Sets Bl and B2). The dotted line is proportional 
to Rm. The inset shows the normalized kinetic helicity of the flow. 



0.02 



0.00 - 



-0.02 



-0.04 



-0.06 




Sh = -2.5 

Sh = -0.25 

Sh = -0.075 

Sh = -0.025 



4- 



0.1 



1.0 10.0 
Rm 



100.0 



10" 
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10" 



10" 
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Rm 



Rm = 0.5 



Rm = 0.2 



Rm = 0.04 



0.01 



■fB" 



0.10 



-Sh 



FIG. 5. Symmetric contribution Q12 as a function of Rm for Pm = 1 
and four values of shear as indicated by the legend (Sets B 1 to B4). 

values of the shear parameter, Sh = —2.5 (—0.25), there is a 
vorticity dynamo for Rm > 1 (Rm > 3), so no points are 
plotted in those cases. 

The off-diagonal component ol\2, shown in Fig.|5] is pro- 
portional to Re for small Rm, while the analytical expres- 
sion ( fTsT i yields 0112 cx Re^. A sign change occurs for 
Rm PS 2, and the values of ol\2 are consistently negative in 
this regime in agreement with Eq. (l25l l derived for Rm ^ 1 
and Re ^ 1. The data is noisy but suggest that a\2 could be 
independent of Rm at high Rm in an agreement with the ana- 
lytical result (I25I 1. Furthermore, for small Rm the dependence 
on shear is weak, although a clearer dependence on shear is 
seen for Rm greater than around 10. 



FIG. 7. Pumping velocity 7 — ^(021 — 012) normalized by ao as a 
function of Sh for Pm = 1 and different values of Rm as indicated 
in the legend (Sets C1-C4). Analytical results according to Eqs. l ll4b 
with Ci(g) = 1, and l l24t are overplotted with dotted and dashed 
lines, respectively. 



In Fig. |6] we show a as a function of Rm. We find that a 
is proportional to Rm for small magnetic Reynolds numbers 
in agreement with Eq. (fTST l. For Rm greater than roughly five, 
a decreases slightly, while the theory suggests that a is inde- 
pendent of Rm for Rm ^ 1. This inconsistency can be un- 
derstood in terms of the relative kinetic helicity 'H.j (kiU^^^, 
where H = a; • tt, which decreases by about 20 per cent be- 
tween Rm 8 and 83 (see the inset in Fig. |6]). Since a (x H, 
this explains the decrease of a with Rm for Rm ^ 1. 
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FIG. 8. Symmetric contribution ai2 normalized by qq as a function 
of Sh for Pm = 1 and different values of Rm as indicated in the 
legend (Sets C1-C4). Runs with Rm = 16 are shown with asterisks 
and connected by a dashed line. For these runs Q12 < so the plot 
shows — ai2. 
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0.010 
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1.000 
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FIG. 9. Pumping coefficient 7 = \{a2\ — 012) as a function of the 
Mach number for Pm = 1 (Set Dl). The normalization factor is 



ao = 



1 , 



rltr 



, and Sh= -0.10. 



3. Dependence on shear 

Figure |7] shows the pumping velocity 7 normalized by ao 
as a function of the shear number, Sh, for Pm — 1 and differ- 
ent values of Rm. Linear dependence of 7 on shear is clearly 
seen in Fig. [T] This is in agreement with the analytical re- 
sult of Eq. ( fT4] l. Rather surprisingly, the data for ai2 suggest 
that there is no dependence on shear (Fig.|8), in contradiction 
with the analytical result of Eq. (flST l that was derived for small 
shear, St^ <^ 1. 

Note that our theory has been developed for incompressible 



flow since the DNS results are nearly independent of Mach 
number for Ma < 0.05. This is shown in Fig.|9] where we no- 
tice a sharp decline of 7 for larger values of the Mach number. 
We are not aware of similar findings for mean-field transport 
coefficients as a function of Mach number. 



IV. DISCUSSION AND CONCLUSIONS 

To clarify the physical effect related to the pumping ve- 
locity, 7, and the off-diagonal components of the tensor Uij 
we rewrite the contributions to the mean electromotive force 
which are proportional to the mean magnetic field in the fol- 
lowing form: 



?{S) 



OitjBj + (7 X B) 
(P) X B^P 



T 



(50) 



where B*^^) is the toroidal mean magnetic field directed along 
the mean shear velocity U (along the y axis), B^^' is the 
poloidal mean magnetic field directed perpendicular to both, 
the mean shear velocity U and the mean vorticity (along the 
X axis), while the pumping velocities, 7^-'") and 7^^', of the 
toroidal and poloidal components of the mean magnetic field 
are given by: 



'■^^ = z {ai2 + 7), 



= -~z {ai2 - 7). 



(51) 
(52) 



Here we take into account the following identities for the off- 
diagonal components of the tensor a,y = {xiVj + Xjiji) ai2 
and oiijBj ~ ai2 z x (B*^^' — i?^"^-*), where ai2 — a2i and 
X, y, z are the unit vectors directed along x, y and z axes, 
respectively. 

It follows from these equations that, when ai2 > 7 > 0, 
the effective pumping velocity of the poloidal mean magnetic 
field is directed upward (along the z axis), while the effec- 
tive pumping velocity of the toroidal mean magnetic field is 
directed downward. When ai2 < 0, but |ai2| > 7, the sit- 
uation is opposite, i.e., the effective pumping velocity of the 
toroidal mean magnetic field is directed upward, while the ef- 
fective pumping velocity of the poloidal mean magnetic field 
is directed downward. Therefore, the effective pumping ve- 
locity, 7, as well as the off-diagonal components of the tensor 
aij, result in a separation of toroidal and poloidal components 
of the mean magnetic field. This effect is very important for 
large-scale dynamo action in shear flow turbulence. 

Another reason for the different pumping velocity of 
toroidal and poloidal components of the mean magnetic field 
is a combination of the effects of rotation and stratification on 
small-scale turbulence. The effect of the separation of toroidal 
and poloidal components of the mean magnetic field was early 
identified in analytic calculations of rotating stratified turbu- 
lence in |26, 56], confirmed in DNS of rotating stratified con- 
vection It57i,i58il . and included in numerical mean-field mod- 
eling of the solar dynamo in lf59ll . Note also that a nonlinear 
feedback of the mean magnetic field to turbulent fluid flow 
causes a different pumping velocity of toroidal and poloidal 
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components of the mean magnetic field ^ . The latter effect 
was included in numerical mean-field modeling of the solar 
dynamo in [60]. 

In summary, we have developed a mean-field theory for a 
pumping effect of the mean magnetic field in homogeneous 
helical turbulence with imposed large-scale shear. In our anal- 
ysis we use the quasi-linear approach, the path-integral tech- 
nique and tau-relaxation approximation, which allow us to de- 
termine all components of the a tensor in different ranges of 
hydrodynamic and magnetic Reynolds numbers. The pump- 
ing effect depends on the a effect and on shear. Using DNS 
and the kinematic test-field method we were able to determine 
all components of the a tensor from numerical simulations 
of sheared helical turbulence. The major part of the numer- 
ical results for the effective pumping velocity, the diagonal 
and off-diagonal components of the a tensor are in a good 
agreement with the theoretical results. However, the numeri- 
cal results for ai2 suggest that there is no dependence of the 
off-diagonal component on shear in contradiction with the an- 
alytical result. In addition, according to the numerical results 
ai2(Re) is proportional to Re for small Rm, while the theory 
yields Q!i2 oc Re^. On the other hand, the change of the sign 
of ai2 from positive for small Rm to negative for large Rm 
observed in DNS is in agreement with the theoretical predic- 
tions. 
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Appendix A: The integrals of the Green functions 

For the integration in cj-space in the case of small magnetic 
and hydrodynamic Reynolds numbers we used the following 
integrals in Eqs. ( fTTT i and (1% : 



loik) = ^ 


f Gf} Gjj G ^ duj — 


h{k) = ^ 


I C* duj ~ 




1 G^G, {Glfduj 


h{k) = ^ 


1 G.,G, {Glfduj 



TT 

TT (5i^ + 77) 



h{k) = 
h{k) = 
Ie{k) = 
Irik) = 
Is(k) = 



J G„ G, (G:)2 doj = 



Gj^ G ^ G ^ did 



TT {3u -h 77) 
TT 



Grj G^ G ^ duj 



G^ Gjy G jy duj 
G Gy G jy duj 



2iy2 [V + T]) fc6' 
TT 

Av^ {v + rj) fcS' 

TT 

v {v + viY k^ ' 

TT 

V (y + 77)2 k^ 



Appendix B: Derivation of Eqs. ( |26t and ^30) in path-integral 
approach 



To derive Eq. ( |26] l we use an exact solution of the induction 
equation with an initial condition B{t = to, x) = B{to, x) 
in the form of the Feynman-Kac formula: 



B,it,x) = {G^,it,to,OBj{to,0)u, , 
and assume that 

-Bj (to,l)= J exp{i^ ■ q)Bj{tQ,q)dq 

Substituting Eq. ( IB2I ) into Eq. ( IB lb we obtain 

Bi{t,x) = J {Gtj{t,to,^) exp[ii-q])y 
X Bj{to, q) exp(iq ■ x) dq , 



(Bl) 



(B2) 



(B3) 



where ^ = ^ — a;. In Eq. ( |B3t we expand the function exp[«^ • 
q] in Taylor series at q = 0: 

00 ^ 



k=0 



and use the identity: 

V'^ exp(ia; • q) ~ (iqf' exp(ia; • q). 
This allows us to rewrite Eq. (IB3l l as follows: 



00 

fe=0 

X y Bj{tQ,q)exp{iq ■ x) dq . (B4) 

After the inverse Fourier transformation, Bj{to,x) = 
J Bj(to,q) exp{iq ■ x)dq, in Eq. (IB4I) we obtain Eq. ( l26l ). 
Equation (IB2I ) can be formally considered as an inverse 
Fourier transformation of the function Bj{to,^). Equa- 
tion ( |26] l has been also derived by a rigorous method, using the 
Feynman-Kac formula and Cameron-Martin-Girsanov theo- 
rem (see il). 

Averaging Eq. (|26] | over the random velocity field yields the 
equation for the mean magnetic field 



B,{{m + l)T,x)^ {{G,,{t,s,^) cxp(^ V))), 
xBj{inT, x), 



(B5) 
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where the angular brackets (•) denote the ensemble average 
over the random velocity field. Now we use the identity 

B,{t + T,x)^exp(^T^^B,{t,x), (B6) 

which follows from the Taylor expansion 



f{t + r) = 
Therefore, Eqs. 



d 



E 

rn—l 



dt 



f{t) = exp 



yield 



+Cij„i„V,„V„)Sj = cxp(tX)B, 



dt I ml 



(B7) 

where G,, = (jG.y))„ = S,, +U,^, t + 0[{St)% = 
iii))^ = ^-U,T + 0[{St)% = ((UG,j))„, 

Cijmn — {{£,m£,nGij)).u,, and wc introduced the operator L, 
which allows us to reduce the integral equation (|B5t to a par- 
tial differential equation. Indeed, Eq. (IB7b . which is rewritten 
in the form 



exp 



B = B, 



reduces to 



dB 

'dt 



LB. 



(B8) 



(B9) 



Taylor expansion of the function exp(TL) reads 

exp(rL) = ^ + tL+ (rL)V2 + (BIO) 
where E is the unit operator. Thus, Eqs. (IB7l i and (IB 10b yield 



L = L, 



-{G. 



I] - Sij 



A 



(BID 

where £)y>„„ = (Gy>„„ - AikmAkjn) /'^t. This yields Eq. 



+ 0(V3), 



Appendix C: Orr- Kelvin random shearing waves for small 
hydrodynamic Reynolds numbers 

We explain here the details that led to the derivation of 
Eqs. ( [36l ) and dJTl ). We seek the solutions of the linearized 
Eq. (|2]l for incompressible velocity field u as superpositions 
of the Orr- Kelvin shearing waves: 

u{t,r) — / u(t,ko)exp[ik{t) ■ r]dko, (CI) 



(see, e.g., 1231 |48!-l50tl), where fco = ikM,ky,k,), k{t) = 
{kxo — Skyt, ky,kz) and we neglected weak Lorentz force. 
The amplitudes of the shearing waves satisfy the following 
equations: 



du,,{t,ko) 
dt 

duz{t,ko) 
dt 



^gkyk^{t) _^,^2 



= 25* 



ky kz 



l^k^t) 



u^it^ko) + fx 



(C2) 



{t,ko)-iyk'^{t) u,{t,ko)+f. 



(C3) 



These equations were obtained by taking twice curl of Eq. (|2|i. 
Equations iC2i and ( |C3b have explicit solutions: 



Ux{t,ko) 



kHt) Jo 



dt'GAt,t')k\t')Ut',ko), 



u,it,ko)^ui'\t,ko) + u<^^Ht,ko), 



(C4) 
(C5) 



j{t, fco) = -1— [kx{t) Ux{t, fco) + ^2 U2(t, fco)] , (C6) 



«(i,fco)= / dt'G,{t,t'),Ut',ko), 



(CI) 



i<^^\t,ko)^2Skyk, dt'^^^^Ux{t',ko), (C8) 



where G,,(t,i') = exp -f dt"k^{t") . Equations (|C4]l- 



for a white-in-time forcing yield the following formulas 
for non-instantaneous two-point correlation functions: 



(^,(t,fco)<«(t',fco)> - G.(t,t') ^ 

X {ux{t',ko)ul'-'\t',ko)), (C9) 
{u<^P{t,ko)K{t',ko)) =G,{t,t') 

X (uW(t',fcoX(i',fco)), (CIO) 

{uxit,ko)u*J^\t',ko))=2Skykz dt"^^0^ 

x{ux{t,ko)ul{t'\ko)), (Cll) 

{ui'Ht,ko)uUt',ko))^2Skykz J^t"^^^ 

x{ux{t",ko)ul{t',ko)), (C12) 



where for t" < t' 



X (u,(i",fco)<(t",fco)), (C13) 



and for t" > t' 



{uxit",ko)uUt',ko)) ^G,it",t')^^ 

X K(i',fco)<(<',fco)). (C14) 
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